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Investigation of the radio wavefront of air showers with 
LOPES measurements and CoREAS simulations 
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Abstract. We investigated the radio wavefront of cosmic-ray air showers with LOPES measurements and CoREAS simula¬ 
tions: the wavefront is of approximately hyperbolic shape and its steepness is sensitive to the shower maximum. For this study 
we used 316 events with an energy above 0.1 EeV and zenith angles below 45° measured by the LOPES experiment. LOPES 
was a digital radio interferometer consisting of up to 30 antennas on an area of approximately 200 mx 200 m at an altitude 
of 110m above sea level. Triggered by KASCADE-Grande, LOPES measured the radio emission between 43 and 74MHz, 
and our analysis might strictly hold only for such conditions. Moreover, we used CoREAS simulations made for each event, 
which show much clearer results than the measurements suffering from high background. A detailed description of our result 
is available in our recent paper published in JCAP09(2014)025 [1], The present proceeding contains a summary and focuses 
on some additional aspects, e.g., the asymmetry of the wavefront: According to the CoREAS simulations the wavefront is 
slightly asymmetric, but on a much weaker level than the lateral distribution of the radio amplitude. 

Keywords: radio detection, cosmic rays, air showers, LOPES 


1 Corresponding author 

Email address: frank.schroeder@kit.edu 



INTRODUCTION 


There are several motivations bringing the wavefront of the radio signal emitted by air showers into the focus of 
interest: First, the wavefront is of principle interest, since the wavefront is linked to the process of the radio emission. 
Thus, understanding the wavefront means understanding the relevant parts of the emission process. Second, the 
wavefront can be used for the reconstruction of the cosmic-ray composition, since the steepness of the wavefront 
is correlated to the distance to the shower maximum, which itself statistically depends on the type of the primary 
particle. Third, the wavefront can be a technical instrument to generally improve the radio detection of air-showers, 
e.g., by a better reconstruction of the shower geometry, or by a better discrimination of air showers against background 
events. 

For the analysis of the radio wavefront we used 316 LOPES events triggered by the co-located KASCADE array 
[2], which fulfill the following criteria of the KASCADE reconstruction: energy E > 10 17 eV, zenith angle 0 < 45°, 
core position within 90m from the center of KASCADE. Thus, for all those selected events there are antennas in 
different azimuthal directions, and the asymmetry of the wavefront partially averages out. Moreover, the events had to 
show a clear radio signal in the cross-correlation beam [3], and we excluded thunderstorm events [4], The analysis of 
the LOPES events followed our standard reconstruction pipeline, with the one exception that we assumed a hyperbolic 
wavefront for the beamforming instead of a spherical or conical wavefront. See references [1, 5, 6] for more details on 
the standard LOPES analysis. 

This means that our results are strictly valid only for LOPES conditions, i.e., in particular for the east-west 
polarization component of the radio signal in the effective frequency range of 43 — 74 MHz, the LOPES altitude of 
110 m above sea level, the geomagnetic field at Karlsruhe, Germany, and the distance range up to 200 m. Nevertheless, 
a preliminary analysis based on a CoREAS simulation indicates that the wavefront has still the same, approximately 
hyperbolic shape when extending the bandwidth from very low frequencies up to 100MHz. 

For each LOPES event we produced two CoREAS simulations [7] (one with proton, one with an iron nucleus as 
primary particle) using the KASCADE reconstruction of the air shower as input parameters. In addition, we simulated 
a few vertical showers with an artificial grid of antennas which are aligned exactly in east, west, north and south 
direction, to study the asymmetry of the wavefront due to the interference of the dominant geomagnetic and the sub¬ 
dominant Askaryan effect. For the CoREAS simulations the geometry is known exactly. Thus, unlike to the LOPES 
measurements, we have not applied the beamforming method to reconstruct the wavefront, but instead performed a 
simple fit to the pulse arrival time at each antenna position (time of the maximum in the east-west component of the 
radio signal filtered to the effective bandwidth of LOPES). 


RESULTS 

The CoREAS simulations show that the radio wavefront is asymmetric in a similar manner to already known 
asymmetry of the lateral distribution [8, 9], but on a much weaker scale as shown in figure 1. Thus, for most purposes 
a symmetric wavefront is a good approximation, especially, since the asymmetry will partially average out in real 
measurements which have antennas in several azimuthal directions from the shower core. 

Therefore, we compared only different models for symmetric wavefronts with each other. In the simulations, the 
hyperbolic wavefront is clearly favored versus the spherical and the conical wavefront, although a conical wavefront 
might be a sufficient approximation for larger distances to the shower axis (d > 50m). The measurements are 
compatible with this result, but the situation is not so clear, since there is significant background. Nevertheless, the 
hyperbolic wavefront is slightly favored by the measurements, too. For example, the average deviation between the 
arrival directions reconstructed by KASCADE and by LOPES is smallest when using the hyperbolic assumption for the 
beamforming. It is (0.683 ±0.026)° when assuming a spherical wavefront, (0.637 ±0.022)° for a conical wavefront, 
and (0.622 ±0.023)° for the hyperbolic wavefront. 

The symmetric, hyperbolic wavefront is described by the following formula for the arrival time T relative to the 
arrival time at the shower core: 


ct(d,z s ) = \J{dsinp) 1 + (c ■ b) 2 + z s cosp + c ■ b 


( 1 ) 


where d is the orthogonal distance to the shower axis, z s is the orthogonal distance to the plane (with z s = 0 at the 
shower core), c is the speed of light, p is the angle between the asymptotic cone of the hyperboloid and the shower 
plane, and b is the offset between the hyperboloid an the apex of the asymptotic cone. Thus, the hyperbolic wavefront 
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FIGURE 1. Asymmetry of the wavefront (left) and the lateral distribution (right) of a vertical shower initiated by a 10 17 eV 
proton simulated with CoREAS in different azimuthal directions from the core, and hyperbolic fits to the wavefronts. 


has two parameters, namely b and p. Both parameters are correlated to some degree, and we fixed b to —3 ns for the 
reconstruction of the wavefront. By fixing b, we improve the correlation of the remaining free parameter p with the 
distance to the shower maximum: The larger the distance to the shower maximum, the flatter the wavefront and the 
smaller p. 

The average values reconstructed for p are: (0.0210 ± 0.0055) rad for the LOPES measurements, (0.0216 ± 
0.0036)rad for the proton simulations, and (0.0194± 0.0030)rad for the iron simulations, i.e., the wavefront devi¬ 
ates from the shower plane by more than 1°. 

The distance to the shower maximum depends on its atmospheric depth X max and the zenith angle 9 of the showers. 
To reconstruct X max from the cone angle p, the zenith angle dependence of p has to be corrected. This correction 
is difficult to calculate, because in first order the shower is not a point source, but an extended line source (which 
only is an approximation, too, because of the lateral extension of the shower). Thus, we have determined the zenith 
dependence of p empirically, and correct for it with a simple parametrization. After the correction, the cone angle p is 
approximately linearly correlated with I max : 

2fmax ~ 25,000g/cm 2 • p/rad • cos~ 1,5 9 (2) 

The coefficients have been determined from the CoREAS simulations: The coefficients for proton and iron simulations 
are similar within a few percent, but not identical. The zenith dependence is the same for LOPES and CoREAS (within 
uncertainties). However, the proportionality coefficient cannot be checked experimentally, since LOPES/KASCADE 
do not feature any reference measurement of A max . 

Figure 2 shows the comparison of the reconstructed and the true A max for the simulations. The reconstruction is 
precise to approximately 25 g/cm 2 . This means that under ideal conditions (no background, known shower geometry), 
light and heavy nuclei can be well separated. When applying the reconstruction derived from the simulations to the 
measurements, the result seems to be reasonable, but the precision of the measurements is significantly worse: most of 
the width in the distribution of the LOPES X max values is due to measurement uncertainties. Since LOPES features a ns- 
precise time calibration [10], we conclude that the main problem is the external radio background at the experimental 
site close to Karlsruhe, which causes significant measurement uncertainties for the pulse time [11], 

The remaining uncertainty of 25 g/cm 2 in the simulations causing the spread in figure 2 (left) is not due to the fixed 
the offset parameter b. But it might be due to the asymmetry of the wavefront, or because even a hyperboloid is not the 
exact shape of the wavefront, but only an approximation. Consequently, a more detailed investigation of the wavefront 
might be necessary when aiming at a A max resolution of better than 25 g/cm 2 . 

CONCLUSION 

As seen clearly for the CoREAS simulations, as statistically indicated for the LOPES measurements, and in agreement 
with recent LOFAR results [12], the radio wavefront of air showers can be better described by a hyperboloid than by 
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FIGURE 2. Left: correlation of the cone angle p after correction for the shower inclination with the true atmospheric depth of 
the shower maximum A max . Right: X max reconstructed by making use of the linear correlation in the left plot (published in [1]). 


a plane, a sphere, or a cone. Although the wavefront is slightly asymmetric, according to the simulations, a symmetric 
wavefront is a sufficient approximation, at least for the currently achieved measurement precision and for distances 
< 200m to the shower axis. 

The steepness of the wavefront, namely the cone angle p can be used to reconstruct A max and to statistically study 
the cosmic-ray composition. The accuracy likely can be improved by combining the wavefront method with other 
independent measurements of the same events. For example, the slope of the lateral distribution of the radio amplitude 
is also sensitive to A max [13, 14]. 

Finally, the improved knowledge of the wavefront leads to a better reconstruction of the shower geometry. Provided 
sufficient measurement precision for the arrival times, wavefront measurements could be used to distinguish air 
showers from other radio sources, e.g., air planes: because many background sources are in good approximation 
point sources, background events usually should have a spherical wavefront and not a hyperbolic one. 
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